PROCEEDINGS 

0F SCIENCE 



Mass Assembly of High-z Galaxies with MASSIV 



V. Perret , B. Epinat, P. Amram, 0. Le Fevre, C. Lopez-Sanjuan, L. Tasca, L. Tresse 

Aix Marseille Universite, CNRS, LAM (Laboratoire d'Astrophysique de Marseille) UMR 7326, 

13388, Marseille, France 

E-mail: valentin . perretSoamp . f r 

^ T. Contini, C. Divoy, J. Queyrel, J. Moultaka 

(N IRAP, Universite de Toulouse, UPS-OMP, CNRS, Toulouse, France 

O 

q D. Vergani 



in 
c3 



(N 



INAFIASFBO, Via P. Gobetti 101, 1-40129 Bologna, Italy 

F. Bournaud 



O CEA, IRFU, SAp, F-91191 Gif-sur-Yvette, France 

^ B. Garilli, L. Paioro 

Q_i INAF-IASFMI, Via E. Bassini 15, 1-20133 Milano, Italy 

M. Kissler-Patig 

ESO, Karl-Schwarzschild-Str. 2, D-85748 Garching b. Munchen, Germany 



MASSIV (Mass Assembly Survey with SINFONI in VVDS) is a sample of 84 distant star-forming 
galaxies observed with the SINFONI Integral Field Unit (IFU) on the VLT. These galaxies are se- 
lected inside a redshift range of 0.8 < z < 1.9, i.e. where they are between 3 and 5 billion years 
£f) old. The sample aims to probe the dynamical and chemical abundances properties of represen- 

tative galaxies of this cosmological era. On the one hand, close environment study shows that 
about a third of the sample is involved in major mergers. On the other hand, kinematical analysis 
£N| revealed that 42% of the sample is rotating disks, in accordance with higher redshift samples. The 

remaining 58% show complex kinematics, suggesting a dynamical support based on dispersion, 
m ^ and about half of these galaxies is involved in major mergers. Spheroids, unrelaxed merger rem- 

nants, or extremely turbulent disks might be an explanation for such a behavior. Furthermore, the 
spatially resolved metallicity analysis reveals positive gradients, adding a piece to the puzzle of 
galaxies evolution scenarios. 
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1. Mass Assembly of High Redshift Galaxies 

In the current framework of ACDM cosmology, the mean energy density of the Universe is 
shared between dark energy, dark matter and baryonic matter. In this context, the formation of 
structures in the Universe is seeded by small perturbations in matter density expanded to cosmo- 
logical scales by inflation. Given such a cosmological model, the Universe is about 13.7 billion 
years old, and galaxies are expected to form around z=20-50 (Gao et al. 2007), when the first suffi- 
ciently deep dark matter potential wells are made and allowed gas to cool and condense to produce 
primeval stars and galaxies. 

This standard model reproduces well the linear initial conditions, the intergalactic medium 
structure during galaxy formation, and large scale structure as observed today. However, the hi- 
erarchical dark matter halo formation paradigm remains strongly debated. The physical processes 
responsible for mass assembly at early epochs is still unclear. While merging events are expected 
to play an important role in the building of the Hubble sequence, smooth cold gas accretion might 
also have strongly contributed to the growth of galaxies. Secular processes such as stellar feed- 
back may also have a major role in the build-up of local galaxies. High-z galaxies tend to have a 
higher gas fraction, thus a powerful stellar feedback would be able to drive the evolution within the 
gravitational potential. 

This debate can be 
addressed with spatially 
resolved measurements 
from 3D spectral anal- 
ysis, allowing to probe 
various physical quanti- 
ties, e.g. star formation, 
gas kinematics, chemical 
abundances. The MAS- 
SIV sample investigates 
a redshift range unex- 
plored by Integral Field 
Spectroscopy (IFS), con- 
sidering that four dis- 
tinct 3D surveys are pre- 
ceding it. IMAGES 
(Puech et al. 2008) sam- 
ple, with 0.4 < z < 
0.75, showed that regular 
rotating disks are quite 
similar to the local ones and that mergers are playing an important role in galaxy mass assembly. 
SINS (Forster Schreiber et al. 2009) sample, Wright et al. (2009) and Law et al. (2009) samples, 
and LSD/ AMAZE (Maiolino et al. 2008) sample, with z>1.5 (mostly with z>2), showed that a lot 
of young galaxies are experiencing a high gaseous turbulence. 

The MASSIV sample was built to probe representative galaxies in the redshift range 0.8 < z 
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Figure 1: Evolution of the cosmic star formation rate density as a 
function of look-back time and redshift. Major IFU surveys redshift 
ranges are compared. The relative height of each box is proportional 
to the sample size. Adapted from Contini et al. (2012). 
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< 1.9 (Fig. 1), where the cosmic star formation rate history is expected to peak (Cucciati et al. 
2012), and where we look forward to a transition between small and disturbed galaxies towards the 
Hubble sequence. This redshift range is also an opportunity to infer the establishment of the stark 
dichotomy within the galaxy population, already in place at z=l. 

MASSIV is an ESO large program (200 hours) 
grouping 84 galaxies observed with the NIR-IFU SIN- 
FONI on the VLT from 2008 until completion in 201 1 
(Contini et al. 2012). The J- or H-bands have been 
used to target the redshifted Ha emission line with a 
high spatial resolution (<0.8"), and a total integration 
time varying between 80 and 120 min. Among these 
84 galaxies, 11 were observed with the adaptive op- 
tics system (AO), reaching a spatial resolutions close 
to 0.20". As of today, it is the largest sample of high-z 
galaxies observed with IFS. 

The strength of this survey lies in its well defined 
parent sample, that is the VVDS (Le Fevre et al. 2005) 
a redshift survey selected in magnitude (Iab < 24) in- 
cluding 35000 spectra in the visible, avoiding any bi- 
ases linked to a priori color selection techniques. A 
high completeness of the parent sample is mandatory if 
one wants to probe normal and representative galaxies. 
Galaxies were selected on [Oil] 3727 equivalent width 
strength or rest-frame UV intensity from SED fitting, 
both being a proxy for star formation. This star forma- 
tion criterion ensures that the brightest rest-frame op- 
tical emission line Ha ([OIIIJ5007 for a few galaxies) 
is available to probe resolved kinematics and chemi- 
cal abundances down to galaxies with a SFR close to 1 
M & yr~ 1 . Galaxies were selected to have a sufficiently 
close bright star valuable for AO/LGS. The continuum I-band magnitude of each galaxy is esti- 
mated using the best-seeing CFHT-LS images with a resolution better than 0.65". 

Fig. 2 compares the relation between stellar mass and star formation rate for the major IFU 
samples, over-plotting empirical relations for different redshifts. On the one hand, we see that for 
a given stellar mass SINS, LSD/AMAZE are globally probing galaxies with higher star formation 
rate than MASSIV, while on the other hand, for a given stellar mass, the IMAGES sample is 
globally probing galaxies with a lower star formation rate than what is expected for a representative 
sample. 

2. Kinematical analysis of MASSIV galaxies 

The full MASSIV sample allows resolved velocity measurement for ~90% of the galaxies. 
Each velocity field was fitted using a PSF-convolved 'flat' model rotation curve (Epinat et al. 2010). 
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Figure 2: SED-derived star formation 
rate as a function of stellar mass. The 
lines represent the empirical relations 
between SFR and stellar mass for dif- 
ferent redshifts between z = and z 
= 3 following the analytical expression 
given in Bouche et al. (2010). All 
the illustrated IFU samples are rescaled 
to match the MASSIV sample redshift 
range. Adapted from Contini et al. 
(2012). 
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A kinematical classification of the first epoch sample has been performed in Epinat et al. (2012), 
using three distinct estimators. A first parameter used to distinguish fast rotators from slow rotators 
is the total velocity shear V s / iear measured on the velocity field, without any inclination correction 
(see Fig. 3). With this simple parameter applied on the full sample, 33 (42%) galaxies exhibit high 
velocity shear (V s /j ear >100 km.s~ l ), and 45 (58%) galaxies have a low velocity shear (V s i iear <l00 
km.s~ l ). 

We also compared a parameter measuring the discrepancy between the position angle of the 
major axis (PA) of the stellar component and the kinematical PA of the gaseous component with 
the mean weighted velocity field residuals normalized by the velocity shear. Similary it gives 33 
(42%) galaxies classified as "rotating", and 45 (58%) galaxies classified as "non-rotating" systems. 

In addition, the immediate en- 
vironment of each galaxy is probed 
to determine whether the galaxy 
lies isolated or not. We search for 
companions in both the Ha and I- 
band image, with a relative veloc- 
ity lower than 1000 km.s~ l and a 
projected distance lower than 50 
h~ l .kpc. Applying this method, 59 
(76%) galaxies are found to be iso- 
lated, and 19 (24%) are not. 

This work on close environ- 
ment has been pushed further away 
by Lopez-Sanjuan et al. (2012). Defining a close pair as a couple of galaxies with a projected radial 
separation lower than 20 h^ 1 .kpc, and a radial velocity difference lower than 500 km.s~ l , we are 
able to recover a major merger rate Rmm x (1 +z) 3 ' 95 in the MASSIV redshift range. 

Although the origin of galaxies with high V shear is clear (disks with ordered rotation), galaxies 
with low V s h ear are more difficult to interpret. Face-on disks, unrelaxed merger remnants, or star- 
forming spheroids could be an explanation to such a behavior. Nevertheless, the orientation of disks 
with spin vectors randomly distributed could only account for 14% of this low-shear population 
considering disks rotating at V mt =200 km.s~ l . The large fraction of interacting, and non-rotating 
galaxies seems to suggest that an important cosmological mass assembly mechanism is at work 
between redshift 1 and 2. 

3. Fundamental relations with MASSIV 

In Vergani et al. (2012) we investigate the fundamental relations using MASSIV data. With 
dynamical arguments, we derive a gas mass that is on average a fraction of ~ 45% of the dynamical 
mass, assuming no central contribution of dark matter, as in Gnerucci et al. (2011), and consistent 
with lower concentration parameter of haloes at high-z (Bullock et al. 2001). This gas content is 
consistent with the fraction derived using the Kennicutt-Schmidt formulation. 

The evolution of the Tully-Fischer Relation (TFR) is expected to be related both to the con- 
version of gas into stars and to the inside-out growth of dark matter halo by accretion. In Vergani 
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Figure 3: (a) High-shear velocity field (VVDS220584167). 
(b) Low-shear velocity field (VVDS020386743). Contours 
are drawing the Ha flux distribution. Adapted from Epinat 
et al. (2012). 
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et al. (2012) the stellar mass TFR shows a negligible, net evolution in the past 8 Gyrs with a large 
scatter that is reduced, but still remarkable, using the S05 index (S05 = yj 0.5 x vj ot + ofi, Fig. 4). 
We interpret this behavior as an evidence of complex physical mechanism(s) at work in our stellar 
mass/luminosity regime and redshift range. We also conclude a marginal evolution in the size - 
stellar mass and size - velocity relations in which disks become evenly smaller with cosmic time 
at fixed stellar mass or velocity, and are less massive at a given velocity than in the local Universe 
in agreement with cosmological hydrodynamical simulations, e.g. Portinari & Sommer-Larsen 
(2007). 
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Figure 4: Left: Stellar mass content of the MASSIV galaxies compared to the dynamical modeling 
Mdyn- The blue/red points corresponds to rotating/non-rotating galaxies according to Epinat et al. 
(2012) criterion. The empty, dashed symbols show galaxies detected with a low signal-to-noise 
ratio (3 < SNR <4.5). The horizontal and vertical lines are the median values of rotators (point- 
dashed blue line) and non-rotators (dashed red line). The relation M^.,, = +M gas is over-plotted 
with the mean value of the gas fraction (45%) for both rotators and non-rotators. Right: Stellar 
mass TFR at z~1.2 (blue line) based on rotating MASSIV galaxies (blue dots), taking into account 
both the ordered and chaotic motions (Kassin et al. 2001, K01). We over-plot the z ~ 1 slope by 
Kassin et al. (2007) (solid black line), and the z ~ 3 slope by Gnerucci et al. (201 1) (long-dashed 
magenta line). The cyan shaded area shows the la error on the zero-point parameter, the cyan 
dotted line and cyan dashed line are the intrinsic and total scatter, respectively. Adapted from 
Vergani et al. (2012). 



4. Abundance gradients with MASSIV 

The first epoch sample (50 galaxies) enabled an abundances analysis in Queyrel et al. (2012). 
The metallicity of galaxies was derived in taking the [NII]/Ha ratio, with the Perez-Montero & 
Contini (2009) calibration. Among these 50 galaxies, 26 metallicity gradients were measured 
inside spatially distinct annular regions defined by the Ha contour map. While half of the sample 
(14 galaxies) is compatible with a zero metallicity gradient, a quarter (7 galaxies) of the restricted 
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sample exhibits positive metallicity gradients. Among these positive gradients, 4 are classified as 
interacting systems, one might be a chain galaxy, and 2 are flagged as isolated. Such features have 
been observed in local interacting galaxies (Werk et al. 2010); it would be an effect of a fresh 
metal-poor gas infall onto the core of the merger remnant. 

The isolated galaxies detected with positive gradients are more puzzling objects. Such systems 
have been observed in Cresci et al. (2009), suggesting that cold gas accretion onto the center of the 
central regions would be a plausible scenario. The radial abundances of the second epoch sample 
will be published soon, and they will extent the statistics of this first analysis. 

5. Simulations 

High redshift disks have been performed in numerical simulations by Bournaud et al. (2011) 
who have shown that the input of large gas fraction (40-60%) in the initial conditions led to chaotic 
and turbulent gaseous disks. Among different question to tackle in the near future, the understand- 
ing of the underlying processes responsible for the non-rotating galaxies in the redshift range 0.8 
< z < 1.9 led us to simulate high redshift disks from idealized mergers. We are sampling different 
orbital parameters and mass ratio through a set of 20 simulations. These simulations will be pro- 
jected on the sky plane at different evolutionary time using numerous projection angles to mimic 
observational configurations. Using the RAMSES code (Teyssier 2002), we are producing simu- 
lations of mergers with a physical resolution reaching 6 pc, coupled with a new implementation 
handling the feedback of OB-type stars emitting energetic ultraviolet photons. With this data we 
aim to disentangle the different scenarios leading to the formation of a population of galaxies with 
a low V s h e ar, as well as the conditions required for the disk survival. 
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